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ABSTRACT. The paper investigates a natural convection flow past 
an impermeable vertical plate embedded in nanofluid saturated 
porous medium. The base fluid is assumed to have variable vis- 
cosity. The self-similar boundary layer equations are solved us- 
ing standard numerical techniques. The validity of the model is 
checked against other previous works and is found to be in good a- 
greement. The results obtained show that the nanoparticle volume 
fraction and the variable viscosity parameter significantly enhance 
the heat transfer properties from the wall to the fluid saturated 
porous medium. 


1. INTRODUCTION 


The study of convective heat transfer in nanofluids flow through 
porous medium is gaining a lot of attention. Nanofluids are solid-liquid 
composite materials consisting of solid nanoparticles or nanofibers with 
sizes typically of 1-100nm suspended in liquid [14]. On the other hand, 
nanofluids are engineered by suspending nanoparticles in traditional 
heat transfer fluids such as water, oil and ethylene glycol. The nanopar- 
ticles used in nanofluids are typically made of metal (Au, Ag, Cu, 
and Fe), metal oxides (CuO, SiO2, A1203, TiO2, ZnO, and Fe304), 
carbides or carbon nanotubes (diamond, graphite, single/multi wall 
carbon nanotubes) [20]. Nanofluids have attracted great interest re- 
cently because of reports of greatly enhanced thermal properties. For 
instance, when less than 1% of the solid particles (copper or carbon 
nanotubes) are added to the base liquids of lower thermal conductivity 
eg., ethylene glycol or oil, solid particles improve the heat conductivity 
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of the liquid by 40% and 150% respectively. Convectional particle- 
liquid suspensions require high concentrations (>10%) of particles to 
achieve such enhancement. [18] observed that dispersion of nanoparti- 
cles in the base fluid increase the thermal conductivity of the fluid and 
therefore, the presence of porous media enhances the effective thermal 
conductivity of the base fluid . On the other hand, the effective thermal 
conductivity produced by a combination of base fluid, nanoparticle and 
porous medium plays a major role in heat transfer enhancement. Mod- 
eling of convection of heat by a nanofluid flow through a porous media 
is a phenomenon of great interest from both theory and application 
points of view. Nanofluids flow in porous media have numerous ther- 
mal engineering applications such as electronic cooling, vehicle cooling, 
transformer cooling and help to improve thermal properties of oil and 
lubricants [2]. 


The theoretical and applied research in flow, heat and mass transfer 
in porous media has received attention during the past three decades 
({19];{8];[7]; [21];[10]). Since the early work of Darcy in the nineteenth 
century, several and intensive investigations have been conducted on 
the laminar flow and heat transfer through porous medium, [12]. Dar- 
cy’s experiment discovered that, in laminar flow through porous medi- 
um, the pressure drop caused by the frictional drag is proportional to 
the velocity at the low Reynolds number range. The Darcy’s experi- 
ment leads to a law called Darcy law which states that, the velocity of 
flow of a liquid through a porous medium due to difference in pressure is 
proportional to the pressure gradient in the direction of flow. The law is 
a phenomenologically derived constitutive equation that describes the 
flow through a porous medium ({13];[3]). Nanofluids are found to ex- 
hibit higher conductive, boiling point and convective heat transfer per- 
formances compared to base fluids [17]. It has showed that an increase 
in the thermal conductivity of a solid-liquid depends on the shape, size 
and thermal properties of the nanoparticles. Nanofluids, when used as 
coolants can provide dramatic improvements in the thermal properties 
at smallest possible concentrations by uniform dispersion and stable 
suspension of nanoparticles in host fluids [2]. By combining nanoflu- 
id with biotechnology components, nanotechnology can have numerous 
potential applications across a wide range of practical applications such 
as agriculture, pharmaceuticals and biological sensors. The potential 
forms of nanomaterials available for use in biotechnological applications 
includes nanowires, nanofibres, nanostructures and nanomachines [5]. 
The commercialization of nanobiotechnological products seems to have 
a potential future and within next a few years many new products of 
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this nature are likely to be used [17]. Since, nanofluid flow and heat 
transfer in the porous medium has many applications and still a new 
area with potential for engineering as mentioned above, then, the study 
has attracted more attention from the research communities to show 
mathematically what really happen in the system of flow of nanofluid. 


Various evidences of nanofluid flow and heat transfer in the porous 

medium have been reported in the literatures by several researchers. 
The first work on convection flow and heat transfer of nanofluids inside 
a 10.66mm-diameter tube by Park and Cho [2] observed a substantial 
increase in the heat transfer coefficient in the tube regime. The bound- 
ary layer flows of nanofluids have been analyzed by [11] and [6] solved 
the non-similar problem of free convection heat transfer from a verti- 
cal plate embedded in a saturated porous medium with an arbitrary 
varying surface temperature or heat flux. [15]investigated the mixed 
convection flow of non-Newtonian fluid from vertical surface saturat- 
ed in a porous medium filled with a nanofluid. [18] investigated the 
natural convection heat transfer of nanofluid along a vertical plate em- 
bedded in porous medium. 
However the impact of variable viscosity of a natural convection flow 
of nanofluid along a vertical impermeable plate embedded in a porous 
medium has received little attention. Hence in the study presented 
here we investigate the effects of variable viscosity of the base fluid 
with nanoparticles, volume fraction of the nanoparticles and the power 
law wall temperature on steady state two-dimensional natural convec- 
tion flow along a vertical plate embedded on the porous medium. The 
problem has an important application to convective flow on the heated 
plate, especially in the compact design devices that cannot be cooled by 
the use of traditional methods rather than using nanoparticle material 
for enhancing thermal conductivity of the base fluids. Furthermore, 
the paper is organised in the following manner. Section two presents 
the problem formulation, Section three is presented the numerical tech- 
nique. Section four will discuss the results and discussion and finally 
conclusion is made in Section five. 
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FIGURE 1. Physical model and coordinate system. 


TABLE 1. Physical properties of base fluid water and 
copper at 20°C (293K), (Anjali,2011) 














p(kg/m*) cp(J/kgK) | k(w/mK) BIK m?/s 
water | 1000.52 4181.8 0.598 | 210*10-* | 0.143*10-6 
copper | 8954 383.1 385 51*107® | 1.12*1071 


























2. PROBLEM FORMULATION 


The current problem considers the steady convection boundary lay- 
er flow of a nanofluid over an impermeable plate embedded in porous 
medium with temperature dependent viscosity. The impermeable plate 
is placed vertically in the porous medium filled with a nanofluid as 
shown in Fig 1. The system is governed by continuity, momentum and 
energy equations in two dimensional laminar flow of saturated nanoflu- 
id. The Darcy-model is used to describe the flow of nanofluid in porous 
medium with small porosity and low flow along a vertical plate. Let the 
convecting nanofluid and the porous medium be at the local thermo- 
dynamic equilibrium with no viscous dissipation effect to the system. 
Also, it is assumed that the physical properties of the base fluid and 
that of the nanoparticles are constant except the density of the fluid 
at the buoyancy term, so that the Boussinesq approximation holds. 
Furthermore, the base (pure) fluid viscosity varies inversely with tem- 
perature and the flow of nanofluid past a vertical plate is assumed to 
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satisfy slip condition. Following the modified equations studied by [4] 
along with the Boussinesq technique and boundary layer approxima- 
tion, then the modified model becomes; 


Ou Ov 
(1) oe 
(2) lot ye a A Ce 
K Ox nf i 
_ Op 
(3) 0= Dy 
OT OT ƏT 
(4) (Pcp) ny E + o) = Rett a 


Noting that v << u in the boundary layer and this shows that there 
are no body forces (including gravity ) in the y-direction and therefore, 
the force balance in the y-direction is zero. Apart form that, the pres- 
sure forces within the boundary layer is established by the hydrostatic 
condition outside the boundary layer and therefore, the unbalance of 
pressure and gravitational forces within the boundary layer due to the 
variation of density of the fluid creates the fluid motion adjacent to 
the vertical plate. Here, u and v are seepage velocity components a- 
long x and y axes. Noting that, x and y are the cartesian coordinates 
measured along the vertical plate and normal to it, respectively. Fur- 
thermore, T is the temperature of the nanofluid within the boundary 
layer and Hng (T), (pcp) np Keff and (Pcp)epp are the dynamic viscosity 
of the nanofluid, heat capacitance, effective thermal conductivity and 
effective heat capacitance of the porous medium, while p is the pres- 
sure and K is the permeability of the porous medium. The subscript 
”nf” defines all the thermophysical quantities of nanofluid in porous 
medium and these quantities are defined as follows-([{18],[9]); 
kp+2k 5 —20(kp—ky) 


pt D tole er)? Pas = P OP) + O — 9) (08) 
(5) keff = are (Pep) ops et (Pcp) ny + (1 — €) (Pcp)m > 
bing (T) = E a (Pep) ng = (1 — P) (Dep) 5 + 9 G) 


knf = 





In eqn.(5), Y, Bp; Pp, (pcp), and kp are the nanoparticle volume frac- 
tion, volumetric thermal expansion, density, heat capacitance and ther- 
mal conductivity of the nanoparticles in the base fluid, respectively, 
while, By, Pf, (Pcp) F and kr are volumetric thermal expansion, density, 
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heat capacitance and thermal conductivity of the base fluid, respec- 
tively. In the model, it is assumed that the vertical plate satisfies the 
power law temperature Tọ and the temperature far away from the plate 
is at the ambient temperature Tæ. The ambient temperature Ta is as- 
sumed to be less than temperature at any arbitrary reference point 
in the boundary layer formed along the plate. Therefore, under these 
situations the boundary conditions associated with this model are as 
follows: 


(6) v=0,T ST +Cr aty =U, z >0 


(7) u=0,T = Taty > co, t >0 


Furthermore, the temperatute of the plate is a power function of 
distance along the plate measured from its leading edge of the plate 
and is defined as; 


(8) CAT + Crò 


where, C > 0 for heated vertical plate and C < 0 for cooled vertical 
plate, whereas A = 0 for isothermal/constant temperature along the 
vertical plate. The pressure in Eqns (2), (3) is eliminated using the 
cross-differentiation approach and therefore the continuity, momentum 
and heat equations are written as; 


ðu ðv 
o oT 
(10) ay (ng (T) u) = (PB) nf aay 


(11) OT OT kept OT 
u Hv = 
Oy  (Pep)ns OY? 





For the present base fluid, following ([1]; [16]), the dynamic viscosity 
of the base fluid is considered in the following form: 


Hoo 1 





(12) Hs (T) 


TEF) TET) 

where uo is a dynamic viscosity of the ambient base fluid; y is a con- 
stant value, it implies that when y = 0, the dynamic viscosity of the 
base fluid flow remains constant throughout the system. On the other 
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hand a and T, are constants, given by; 


(13) a=, T, — To = —- 
Loo Y 
We now look for similarity solutions of Equations (10) and (11) with 
the boundary conditions (6) and (7) of the following form; 
(14) 
= 1/2 EA (ean i a E 

Y = aes Ra,’ f(n), n = 7 Ras“, O(n) = Tee E 
provided that Ra, = g8 K (Ty, — Tæ )£/Væ@eşp is the local Rayleigh num- 
ber of the porous medium, v, aerr and 0, are the kinetic viscosity, 
effective thermal diffusivity of the porous medium with nanofluid and 
viscosity-variable parameter. Also, WV is the stream function which is 
expressed in terms of derivatives of the seepage velocity as u = OW /Oy 
and v = —OWV/0z, apart from that, the seepage velocities usually satsi- 
fy the continuity equation of the flow. Using, similarity solutions, the 
seepage velocities may be expressed in the following form; 


V 


(15) f 








Gey gas. fos (Ty Te 


yy 

1+A)f-(1-A)nf 
; eee) atas- = at 
where, primes represent the derivative with respect to 7. 
Substituting the seepage velocities Egns (15) and (16) in Eqns (10) 
and (11) together with Eqns (12-14), we get the following system of 
non-linear ordinary differential equations: 


£0 as ( oth fy 9\_ 
an ripa- (etting) (1-7) =0 








/ 


(18) 0 + fo —\fd=0 





»  (1+A) 
2 


Subjected to the boundary conditions: 
y=0, 7 =90, f(0) =90, (0) =1, 


(19) ; 
y > œ, n> œ, f (co) = 0, A(co) = 0 


It has showned that 6,. is determined by the viscosity /temperature char- 
acteristics of the base fluid and the operating temperature difference is 
AT = Ty — T. It is also shown that, if 0, is small then either the base 
fluid viscosity changes with temperature difference or the temperature 
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difference is high. It may be noted that 0, is negative for liquids and 
positive for gases [1]. 


The quantity of physical interest is local Nusselt number Nuz, which 
is defined as follows: 


hx qd. n” OT 
Nuz; = ; h = = ; qu = —ke (=) 
keff Tw =. Too n Oy y=0 


1/2 
dw = hess (Z) = kes0’ (=) a 
y=0 





(20) 


Oy UQef f 


Therefore, using the similarity solutions Egn (14), we get the local 
Nusselt number equation in dimensionless form; 


gBK (To = L 


VQef f 


Nuz , 

(22) aia = [0] 

where by q,, and h are the local heat flux and heat transfer coefficient 
from the vertical plate, respectively. Hence, the effects of the nanopar- 
ticle volume fraction (y), viscosity-variable parameter (6,) and param- 
eter (A) on the velocity profile, temperature profile and convective heat 
transfer (local Nusselt number) from a vertical plate to the nanofluid 
are discussed in the results and discussion part. 





Gi: Nuss | z] i 6'(0) = Rat”? |-6'(0)| 


3. NUMERICAL TECHNIQUE 


The dimensionless flow Equations (17) and (18) are non-linear cou- 
pled ordinary differential equations, which have no exact solution/closed- 
form solution. Thus, they must be computed numerically together 
with the boundary values provided in Egn (19). The solutions of Eqns 
(17) and (18) subject to the boundary condition (19) were obtained 
by applying the fourth-order Runge-Kutta integration with a shoot- 
ing method. In general, the shooting method transforms the original 
boundary value problem into initial value problem. It is an iterative 
technique which attempts to find an appropriate initial conditions of 
the original boundary value problem . Furthermore, it attempts to 
assume some of missing initial conditions which are not provided in 
the original boundary value problem. Therefore, the systems of the 
first order differential equations formed after transforming the orig- 
inal boundary value problem using shooting method were solved by 
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fourth order Runge-Kutta method. The accuracy of the assumed ini- 
tial condition is then checked by comparing the approximated value of 
the dependent variable at the end point of the domain with its known 
value at that end point. If it is happen that the difference still exists, 
then an improvement is done by choosing another initial conditions and 
this process repeated until the difference become very small between 
the calculated and the know value at the end point of the domain. 
Furthermore, the different step sizes were selected to satisfy the con- 
vergence criterion of 1076. From the process of numerical computation, 
the velocity, temperature profiles and local Nusselt number were ob- 
served and their numerical values were presented in the graphical and 
tabular forms. 


TABLE 2. Numerical values of [—6'(0)] for Cu nanopar- 
ticles when y —> 0 and 6, — —oco 



































| [—6"(0)|- Present results | [—6"(0)|- Nield and Bejan (2006) results 
“173 0.0247 0 
-1/4 0.1755 0.162 
0 0.448 0.444 
1/4 0.6286 0.630 
1/3 0.6793 0.678 
1/2 0.7715 0.761 
3/4 0.8931 0.892 
1 1.0005 1.001 

















4. RESULTS AND DISCUSSION 


In our present study, we have computed and analyzed a theoreti- 
cal solution of the dimensionless equations subject to the boundary 
conditions. The computations have been done for the vertical plate 
embedded in saturated porous medium with nanofluid. The tempera- 
ture of the vertical plate assumed to be a power function of distance 
along the plate from its leading edge. Also, the values of the physical 
properties of the base fluid water and copper were taken from ((2]), 
whereas the nanofluid correlation formulas were taken from [9]. The 
effects of nanoparticle volume fraction vy, variable-viscosity paramete 0, 
and parameter À on the dimensionless velocity, temperature and local 
Nusselt number were discussed and shown graphically. Furthermore, 
the present results have been compared with [4] and it is found that 
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they are in good agreement as the variable viscosity parameter goes to 
infinity and nanoparticles in the base fluid approaches to zero, Table 2 
and Fig 9. 


The velocity profile distribution for different values of nanoparticle 
volume fraction y, with fixed parameter \ = —1/3 and variable vis- 
cosity parameter 0, = —1 is depicted in Fig 2. It is observed that 
as the nanoparticle volume fraction increases, dimensionless velocity 
of the nanofluid along a vertical plate decreases. We suggested that 
the decline of the dimensionless velocity of the nanofluid is due to the 
increase of surface drag caused by more addition of nanoparticle in the 
base fluid water. 


In Fig 3, the dimensionless temperature plotted against the simi- 
larity variable 7 for different values of nanoparticle volume fraction. 
From the figure, it is observed that with the increase in nanoparti- 
cle volume fraction, the dimensionless temperature profile increases in 
thermal boundary layer. Here, the nanoparticle volume fraction has 
shown a significant in enhancing heat transfer from the heated plate to 
the fluid in motion adjacent to the plate. 


Figures 4-5 show the dimensionless velocity and temperature profiles 
for different parameter A with fixed nanoparticle volume fraction and 
variable viscosity parameter. It is observed that both of them show a 
common behavior. It means that, as the parameter A increases, veloc- 
ity and temperature profiles decrease. Also, through different compu- 
tations we have observed that the problem becomes realistic when the 
parameter À range is —1/3 <A < 1. 


The dimensionless velocity profile plotted against the similarity vari- 
able 7 was shown in Fig 6, with fixed nanoparticle volume fraction » 
and parameter A. It is observed that as the variable viscosity parame- 
ter 0, decreases, velocity profile increases. This is due to the fact that, 
as the difference in temperature between a vertical plate and ambient 
become smaller, the variable viscosity parameter becomes large, this 
leads to lower the velocity of the fluid. However, when the difference 
in temperature increases the variable viscosity parameter decrease as 
a result the velocity profile of the nanofluid do increase. Furthermore, 
the velocity profiles drop faster as you move away from the vertical 
plate as the variable viscosity parameters decrease. 
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In Fig 7, the effects of variable viscosity parameter 0, on temper- 
ature profiles are depicted. It is observed that, as variable viscosity 
parameter decreases, temperature profiles decrease. The decline of the 
temperature profile is due to the overhauling of heat in the boundary 
layer caused by convection of fluid at high velocity adjacent to the 
heated plate. 


Fig 8 displays the effects of the ratio of local Nusselt number and 
Rayleigh number on the surface of heated plate for various values of 
nanoparticle fraction volume with fixed A and variable viscosity pa- 
rameter 6,. It is observed that, the ratio of local Nusselt number and 
Rayleigh number with nanoparticles are linearly relationship. This in- 
dicates that as nanoparticle increases in the base fluid water, the ratio 
number increases. This indicates that, more heat is removed from a ver- 
tical plate and this action proves that nanofluids are efficient coolants 
than the traditional fluids. 





























FIGURE 2. Effects of nanoparticle fraction volume on 
velocity profiles f(n) when \ = —1/3,0, = —1. 
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FIGURE 3. Effects of nanoparticle fraction volume on 
temperature profiles 6(7) when A = —1/3,6, = —1. 














FIGURE 4. Effects of parameter À on velocity profiles f 
when y = 0.02, 6, = —1. 
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FIGURE 5. Effects of parameter A on temperature pro- 
files 6(7) when y = 0.02, 6, = —1. 





9, = -1000 
9, = -100 
@=-5 
@=-1 








FIGURE 6. Effects of variable -viscosity parameter on 
velocity profiles f’ when \ = —1/3, œ = 0.02. 
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FIGURE 7. Effects of variable -viscosity parameter on 
temperature profiles 6(7) when A = —1/3, @ = 0.02. 











FIGURE 8. Effects of heat transfer from the heated sur- 


face for various values of nanoparticle fraction volume 
when à = —1/3 and 6, = —1. 
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5. CONCLUSION 


Natural convection flow past an impermeable vertical plate embed- 
ded in nanofluid saturated porous medium with temperature dependent 
viscosity is studied numerically using fourth order Runge-Kutta inte- 
gration with shooting method. The numerical solutions of the model 
are compared with previous works and found to be in good agreement. 
The effects of the nanoparticle volume fraction, variable-viscosity pa- 
rameter and the power parameter of temperature along the plate are 
discussed and presented in graphical forms. In general a copper-water 
nanofluid has shown an improvement in heat transfer from a heated 
vertical plate into convective fluid in porous medium. 








T T 
Present results 
Nield and Bejan results 














1.25 














FIGURE 9. Comparison between the present results and 
Nield and Bejan results on the heat transfer from the 
heated plate when ¢ = 0 and 6, + —oo. 
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